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. . Abstract 

Q ' A most important aspect in the field of traffic modeling is the simulation 
of bottleneck situations. For their realistic description a macroscopic multi- 

' lane model for uni-directional freeways including acceleration, deceleration, 

(N 

' velocity fluctuations, overtaking and lane-changing maneuvers is systemat- 

■ ically deduced from a gas-kinetic (Boltzmann-like) approach. The resulting 

> 



equations contain corrections with respect to previous models. For efficient 
computer simulations, a reduced model delineating the coarse-grained tem- 



■ poral behavior is derived and applied to bottleneck situations. 
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I. INTRODUCTION 



Apart from microscopic traffic models, in the last decades a number of interrelated macro- 
scopic traffic models have been proposed |[l|-[TO|. The motivations for developing these were 



to describe and understand the instabilities of traffic flow p|-p|,p!0|-p!2 



to optimize traffic flow by means of on-line speed-control systems JTSj-P^ 



to make short-term forecasts of traffic volumes for re-routing measures P"B|-P^], 

to calculate the average travel times, fuel consumption, and vehicle emissions in de- 
pendence of traffic volume |3,p!9|. 



to predict the effects of additional roads or lanes [p]9|-|22 



Most of these models are restricted to uni-directional freeway traffic and treat the different 
lanes of a road in an overall manner, i.e. like one lane with higher capacity and possibilities 
for overtaking. However, this kind of simplification is clearly not applicable if there is a 
disequilibrium between neighboring lanes. Therefore, some researchers carried out empirical 
investigations of the observed density oscillations between neighboring lanes or proposed 
models for their mutual influences p3|-[28[1. 

However, these are phenomenological models which treat inter-lane interactions in a 
rather heuristic way. Moreover, most of them base on the simple traffic flow model of Lighthill 
and Whitham which assumes average velocity on each lane to be in equilibrium with density. 
This assumption is not very well justified, especially for unstable traffic which is character- 
ized by evolving stop-and-go waves |]5|-p|,p!0HT^ . It is also questionable for lane mergings 
or on-ramp traffic where frequently a disequilibrium occurs PJ^. However, instabilities or 
disequilibria may decrease the freeway capacity considerably. 

Another approach including a phenomenological velocity equation has been proposed by 
Michalopoulos et al. |^. It bases on Payne's model which has been severely criticized 
for several reasons P, p^ , p9| -p^ . Therefore, we will derive a consistent macroscopic multi-lane 



model from a gas-kinetic level of description. This is related to Paveri-Fontana's approach 



(cf. Sec. II), but explicitly takes into account overtaking and lane-changing maneuvers. The 
corresponding Boltzmann-like model allows to deduce macroscopic traffic equations not only 
for the vehicle densities on the different lanes, but also for the associated average velocities 
(cf. Sec. III). Due to different legal regulations, the traffic dynamics on American freeways is 
different from that on European ones (which will be called "autobahns" in accordance with 
Kiihne et al. Hi). 

For efficient computer simulations of large parts of a freeway system it is desirable to have 
a somewhat simpler model. Therefore, in Section IV we will eliminate the velocity equations 
and derive a reduced multi-lane model for the traffic dynamics on a slow time scale. By means 
of computational results it is demonstrated that even the difficult bottleneck situations can 
be successfully simulated with this model. 

A summary and outlook is presented in Sec. V. However, before the macroscopic multi- 
lane model sketched in Ref. is founded, derived, simplified, and simulated, the alternative 



microscopic approaches shall be mentioned. One class of microsimulation models bases on 
cellular automata, like the ones by Rickert et al. [BO] and Nagatani [P7|. These update the 



vehicle dynamics within two successive steps: either the vehicle motion in one step and 
lane-changing in the next step [^, or the left lane in the ffist step and the right lane in 
the second one Bottlenecks were represented by crashed cars with zero velocity |]37 



Noteworthy are also the event-oriented model by Wiedemann and Benz and the social 



force model by Helbing and Schwarz [19 



II. BOLTZMANN-LIKE MULTI-LANE THEORY 

The first Boltzmann-like (gas-kinetic) model was proposed by Prigogine and co-workers 
T|. However, Paveri- Font ana has pointed out that this model has some peculiar 



properties. For this reason, Paveri-Fontana proposed an improved model that overcomes 
most of the short-comings of Prigogine's approach. Nevertheless, his model still treats the 
lanes of a multi-lane road in an overall manner. Therefore, an extended Paveri- Font ana-like 
model will now be constructed. 

Let us assume that the motion of an individual vehicle a can be described by several 



variables like its lane ia{t), its place ra{t), its actual velocity Va{t), and its desired velocity 
foa(t) in dependence of time t. The phase-space density pi{r,v,VQ,t) is then determined by 
the mean number Ani{r,v,vo,t) of vehicles on lane i that are at a place between r — Ar/2 
and r + Ar/2, driving with a velocity between v — Av/2 and v + Av/2, and having a desired 
velocity between vq — Avq/2 and + Avq/2 at time t: 
Ani{r,v,Vo,t) 



pi{r,v,vo,t) = 
1 



Ar Av Avo 

r+Ar/2 v+Av/2 s;o+Ado/2 

T.^nc.{t) fdr'6{r'-r^{t)) fdv'6{v'-v^{t)) f dv',5{v'^ - vl{^)) . 



Ar Av Avn 

r-Ar/2 v-Av/2 vo-Avo/2 



(1) 



Here, Ar, Av, and Awq are small intervals. 5ij denotes the Kronecker symbol and 5{x — y) 
Dirac's delta function. The notation indicates that a time-dependence only occurs in 
exceptional cases. Lane numbers i are counted in increasing order from the right-most to the 
left-most lane, but in Great Britain and Australia the other way round. (For Great Britain 
and Australia "left" and "right" must always be interchanged.) 

Now we utilize the fact that, due to the conservation of the number of vehicles, the 



phase-space density pi{r,v,vo,t) on lane i obeys the continuity equation ||19| , |42|J43 



dpi d 9 (dpi\ (dpi\ (dpA 

+ i>ti.r, V, vo, t) - i>r(r, V, vo, t) . (2) 

The second and third term describe temporal changes of the phase-space density pi{r, v,Vo,t) 
due to changes dr/dt = v of place r and due to acceleration /°, respectively. We will assume 
that the vehicles accelerate to their desired velocity vq with a certain, density-dependent 
relaxation time Ti, so that we have the acceleration law 

mr,v,vo,t) = '^. (3) 

The terms on the right-hand side of equation (^) reflect changes of phase-space density 
Pi{r,v,Vo,t) due to discontinuous changes of desired velocity Vq, actual velocity v, or lane i. 
z/j^(r, V, Vq, t) and i^j/ir, v, Vq, t) are the rates of vehicles entering and leaving the road at place 
r. They are only different from zero for merging lanes at entrances and exits respectively. 



The term 



where 



Pi{r, v,t) = J dvo pi{r, v, Vq, t) (5) 

is a reduced phase-space density and Tj. ^ 1 s is about the reaction time, describes an adapta- 
tion of the actual distribution of desired velocities Poi{vo]r,t) to the reasonable distribution 
of desired velocities pQi{vo; r, /) without any related change of actual velocity v. 

For the reasonable distribution of desired velocities we will assume the functional depen- 
dence 

P,,{vo; r,/) = 1 e-t''"-^-]'/^^^] (6) 



which corresponds to a normal distribution and is empirically well justified [^p4| , p!9| . The 
mean value Voi = ^oil'^)^ and variance 6oi = 6Qi{r,f) of Poi{vo; r,f) depend on road conditions 
and speed limits. Since European autobahns usually do not have speed limits (at least in 
Germany), 6'oj is larger for these than for American freeways. In addition, on European 
autobahns Voj increases with increasing lane number i since overtaking is only allowed on 
the left-hand lane. 

Before we specify the Boltzmann-like interaction term {dpi/dt)i^t and the lane-changing 
term {dpi/dt)\c we will discuss some preliminaries. For reasons of simplicity we will only treat 
vehicle interactions within the same lane as direct pair interactions, i.e. in a Boltzmann-like 
manner 



p5| . Lane-changing maneuvers of impeded vehicles that want to escape a queue 
(i.e. leave and overtake it) may depend on interactions of up to six vehicles (the envisaged 
vehicle, the vehicle directly in front of it, and up to two vehicles on both neighboring lanes 
which may prevent overtaking if they are too close). Therefore, we will treat lane-changing 
maneuvers in an overall manner by specifying overtaking probabilities and waiting times of 



lane-changing maneuvers (which corresponds to a mean- field approach, cf. Ref. [^]). These 
probabilities and waiting times dependent on the vehicle densities and maybe also on other 
quantities. 



For not too large vehicle densities the Boltzmann-like interaction term can be written in 
the form ]19[ 

= YlJ^'"' J ^'^ J dwoW2iv,i\v',i';w,i')pi'{r,v',Vo,t)pi^{r,w,Wo,t) (7a) 

' w<v' 

~Y1 j dv' J dw J dwoW2{v\i'\v,i;w,i)pi{r,v,vo,t)pi{r,w,wo,t) . (7b) 

Term (^a]) describes an increase of phase-space density pi{r, v,vo,t) by interactions of a vehi- 
cle with actual velocity v' and desired velocity vq on line i' with a slower vehicle with actual 
velocity w < v' and desired velocity wq causing the former vehicle to change its velocity to 

V 7^ v' or its lane to i ^ i'. The frequency of such interactions is proportional to the phase- 
space density pif{r,w,Wo,t) of hindering vehicles and the phase-space density pj/(r, t>', t>o, t) 
of vehicles which can be affected by slower vehicles. Analogously, term ( [7E| ) describes a de- 
crease of phase-space density pi{r,v,vo,t) by interactions of a vehicle with actual velocity 

V and desired velocity Vq on line i with a slower vehicle with actual velocity w < v and 
desired velocity wq causing the former vehicle to change its velocity to v' ^ v or its lane 
to i' 7^ i. Since the interaction is assumed not to influence the desired velocities vq, wq, the 
interaction rate W2 is independent of these. However, the interaction rate W2{v' ,i'\v,i;w,i) 
is proportional to the relative velocity |f — w| of approaching vehicles. Therefore, we have 



the following relation: 

W2{v',i'\v,i;w,i) = pf\v - w\6i>(^i+i)6{v' - v) (8a) 

+ p:r\v ~ w\6i>(^i^i^6{v' - v) (8b) 

+ {1 - Pi)\v - w\6iH6{v' - w) . (8c) 



Term (^) describes an undelayed overtaking on lane i' = i + 1 without any change of velocity 
(f ' = v) by vehicles which would be hindered by slower vehicles on lane i. pf denotes the cor- 
responding probability of immediate overtaking. Analogously, term (^) reflects undelayed 
overtaking maneuvers on lane i' = i — 1 with probability p~. Term (|8^) with 

Pi = pf + Pi (9) 
delineates situations where a vehicle cannot be immediately overtaken by a faster vehicle so 




that the latter must stay on the same lane {i' = i) and decelerate to the velocity v' = w oi 
the hindering vehicle. 

We come now to the specification of the lane-changing term {dpi/dt)\c. This has the form 
of a master equation: 

(^) = E W,W)p,{r,v,v,,t) (10a) 

\ / Ic i'i^i) 

- W^{i'\i)p,{r,v,VQ,t). (10b) 

Term ( |10a| ) describes an increase of phase-space density pi(r, t>, t>o, t) due to changes from lane 
i' ^ i to lane i by vehicles with actual velocity v and desired velocity vq. The frequency of 
lane-changing maneuvers is proportional to the phase-space density Pi'{r, v, vq, t) of vehicles 
which may be interested in lane- changing. Analogously, term (|10b| ) refiects changes from lane 
i to another lane i' causing a decrease of pi{r, v, vq, t). For the corresponding rate iyi(2'|2) of 
lane-changing maneuvers we have the relation 

Wi{i'\i) = i^i'Ci+i) + 7^5i'(i-i) , (11) 

since vehicles can only change to the neighboring lanes i' = i ± 1. [Tj~] denotes the 
waiting times for delayed overtaking or spontaneous lane-changing maneuvers on the left- 
hand [right-hand] lane. 

Due to different legal regulations, the explicit form of the overtaking probabilities pf 
and the waiting times in dependence of the vehicle densities is different in American 
countries compared to European ones. A more detailled discussion of this aspect is presented 
in Ref . IB 

III. DERIVATION OF MACROSCOPIC TRAFFIC EQUATIONS 

The gas-kinetic traffic equations are not very suitable for computer simulations since they 
contain too many variables. Moreover, the phase-space densities are very small quantities 
and, therefore, subject to considerable fiuctuations so that a comparison with empirical data 
is difficult. However, the special value of gas-kinetic traffic equations is that they allow a 



systematic derivation of dynamic equations for the macroscopic (collective) quantities one 
is mainly interested in. 

A. Definition of Variables 

The most relevant macroscopic quantities are the vehicle densities 

Pi{r,t) ^ j dv J dvoPi{r,v,vo,t) (12) 

and the average velocities 

Vi{r,t) = {v)i = J dvvPi{v;r,t) (13) 

on lanes i. Here, we have applied the notation 

F,{r,t) ^ {f{v,Vo)h = fdv [ dvof{v,Vo /'^'^''^'Y^ (14) 

J J Pi{r,t) 

and introduced the distribution of actual velocities 

Pi{v] r,t)= / dvo — T — = —r- 15 

J pi{r,t) pi{r,t) 

on lane i. Analogous quantities can be defined for vehicles entering and leaving the road at 
entrances and exits, respectively 

.t(r.t)-JdvJdv,.H':v,v„,t) (16) 

are the rates of entering and leaving vehicles, and 

K^(r,t) = {v)f = I dvvPt{v;r,t) (17) 

their average velocities, where 

Ptiv;r,t)= fdvo%^ (18) 

are the velocity distributions of entering and leaving vehicles, respectively. In addition, we 
will need the velocity variance 

9i{r, t) = {[v - Vi{r, t)f)i ^ j dv[v - Vi{r, t)fPi{v; r, t) = {v^)^ - {{v)if (19) 



and the average desired velocity 



Vo.{r,t) = jdv I dv^v /'^'^'^'Y^ (20) 

Piyi t) 



on each lane i as well as the average interaction rate 
1 1 



of a vehicle on lane i with other vehicles on the same lane. 



dvpi{r,v,t) J dw {v — w)pi{r,w,t) (21) 



B. Derivation of Moment Equations 

We are now ready for deriving the desired macroscopic traffic equations from the gas- 
kinetic equation (0) with (|^), (0), (|^), (|]), (|10]), and (|l5). Integration with respect to vq 
gives us the reduced gas-kinetic traffic equation 

= -(1 - Pi)piir, V, t) J dw{v - w)pi{r, w, t) (22b) 
+ pt_iPi-i{r,v,t) dw {v - w)pi-i{r,w,t) (22c) 



+ Pi+iPi+i{r,v,t) J dw (v - w)pi+i{r,w,t) (22d) 

- {vt + Pi)Piir, V, t) j dw{v- w)pi{r, w, t) (22e) 

1 . 1 . 

+ ^Pi-i(r, w, t) - ^Pi(r, t) (22f) 

+ -^p,+i{r,v,t) - ^pi{r,v,t) (22g) 

+ z>+(r,t;,t)-i>r(r,t;,t) (22h) 



with 



Vo^iv) ^ Vo.iv;r,t) = / dvoVo^^^^f^^ (23) 



and 



(r,t;,t) = / (it;o^>i (^^,^0,^)- (24) 



In formula (^2]), the deceleration term (|22b|) stems from (|8^), the terms ( |22c|) to ( |22e|) 



reflecting immeditate overtaking come from (Ral) and (|8b[), and the lane-changing terms 



( P2i| ), ( |22g| ) originate from (|TT]). The adaptation term {dpi/dt)^^ yields no contribution. 

We will now derive equations for the moments {v^) by multiplying p^ ) with and 
integrating with respect to v. Due to 

= -^((^^'-%).-(A) (25) 

and 

[dvp,ir,v,t) I dw{wv^-v^+^)p,{r,w,t) = {l-p,){p:)\{v),{v'')i-{v^+^),) (26) 



we obtain the macroscopic moment equations 











^i-^Pi{v'')i 

i 


+ + Pi-l{v'')i- 


^i-^^Pi{v^)i 

i 

-1 - ^Pi{v^)i 
^ i 




-1 - -f=Pi{v^)i 

i 


+ z/+(r,t)(t;V- 





(27) 



{v')t = fdv /rf.o.^fi^^^ = fdvv^'-^^ (28) 



Here, we have introduced the notation 

~z/±(r,t) y^^^ z/±(r,t) 

and applied the approximation 

dv pj{r,v,t) / dwv^{v — w)pj{r,w,t) 



w<v 



{v'), I dvp,{r,v,t) I dw{v-w)p,{r,w,t) = (A'^^ (29) 



3 



which is empirically justified due to the smallness of the velocity distributions Pj{v\ r, t) (i.e. 
due to Je'j<^Vj MM)- 



C. Fluid-Dynamic Multi-Lane Traffic Equations 



In order to derive dynamic equations for the densities pi and average velocities V^, we 
need the relations 

{v^)^ = m + {v- Vi)f), = + 2V,{v - V>), + iiv - y,)^), = {y,f + (30) 

and 

Applying these and and using the abbreviations 

' i ^ i -'- i 

equation ( P7| ) gives us the density equations 

^ + ^i% = -P-^ + ''t{r,t)-^i{r,t) (33a) 
+ ^-4 + ^-^. (33b) 

'^i-l '^i ''1+1 '^i 

This result is similar to previous multi-lane models. However, by some lengthy but straight- 
forward calculations we additionally obtain the velocity equations 

dVi dVi dVi Pi 

+ §r^(^.~i - Vi) + §±i(l^,+i - (34b) 



+ '^tiVi'-V^-u7{Vr-V^ (34c) 
with the so-called traffic pressures [^, ^ , ^ 

Vi = pidi (35) 

and the equilibrium velocities 

Vt = V^,-n{l-pi)pA. (36) 



Equation ( p^ ) corrects the phenomenological approach by Michalopoulos et al. The 
terms containing the rates vf and v~ reflect entering and leaving vehicles, respectively. 



Whereas the terms ( |33a|) and ( P4a| ) correspond to the effects of vehicle motion, of accelera- 
tion towards the drivers' desired velocities, and of deceleration due to interactions, the terms 
( P3 b| ) and ( p4b| ) arise from overtaking and lane-changing maneuvers. ( |34b| ) comes from dif- 
ferences between the average velocities on neighboring lanes and tends to reduce them. The 
term ( p4cD has a similar form and interpretation like ( |34bD . It is only negligible if entering 
vehicles are able to adapt to the velocities on the merging lane and exiting vehicles initially 
have an average velocity similar to that on the lane which they are leaving so that V-^ ~ Vi. 

In order to close equations (^) and (|3^) , we must specify the interaction rates 1/T° and 
the variances 6i. Utilizing that the empirical velocity distributions Pi{v]r,t) are approxi- 

we have 



mately normally distributed |10,2 



P,{v;r,t) 



1 ^-\v-V,(r,t)]y\2edr,t)] 



y/2T7ei{r,t) 



(37) 



which implies 



(38) 



With a detailled theoretical and empirical analysis it can be shown P|TP[] that the variances 
6i{r,t) can be well approximated by equilibrium relations Of{pi) which are given by the 
implicit equation 
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(39) 



Here, denotes the overtaking probability for vehicles on lane i, when the densities pj on 
the different lanes j are in equlibrium. For the average desired velocities Voj we have 



Voi = Vo^{r,t)^Voi{rJ) 



(40) 



smce 



Poi{vo] r,/) - PoiivQ] r, t) ^ 



(41) 



due to the smallness of T.. 
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IV. DERIVATION AND SIMULATION OF A REDUCED MULTI-LANE MODEL 



The velocity equations are mainly needed to model the observed traffic instabilities which 
lead to the spontaneous formation of stop-and-go waves at medium densities P-|8|.[T0|-[12 



However, if one is not interested in the density oscillations but only in the average temporal 
evolution of traffic flow, the velocity equations can be eliminated. In order to do this, we will 
apply a method that has been suggested by Sela and Goldhirsch PSI: First, we introduce 

the time averages 

t+AT/2 

Hr,t) = ^ J dtF,{r,t) (42) 

t-AT/2 

over the least common multiple AT of the occuring oscillation periods ATj. Then, the 
quantities p(r, t) and V{r, t) will describe the coarse-grained tiaSic dynamics, in other words: 
the traffic dynamics on a slow time scale. Additionally, the time averages of the total time 
derivatives dpi/dt ad dVi/dt will approximately vanish: 



dPi dpi ^,dpi „ dVi dVi ^^dVi „ 

This corresponds to the assumption that, in coordinate systems moving with velocities 
Vi{r,t), the densities Pi{r,t) and velocities Vi{r,t) oscillate around their (slowly changing) 
equilibrium values. 



Now, we approximate time averages Fi{pi,Vi) of density- and velocity- dependent func- 
tions Fi{pi, Vi) by a series in spatial derivatives of JJl and Vi. For our purposes it is sufficient 



to truncate the expansion after the ffist order |T9|] : 



dpi , ^ ,_ dVi 



Fi{pi, ^ Foo(p-, Vi) + F,o{pl, V,)-^ + Foi(p-, V,)-^ . (44) 
With this and (43) we obtain from the time average of velocity equations 



— \Pi) H J^Vi-l H —Vi+i — 



'~ Pi V£i Pm ■ ^ ' 



"^i ''1+1 



Here, we have restricted our considerations to the case of a freeway without entrances and 
exits. Resolving ( ^5]) with respect to Vi leads to a relation of the form 



V^ = V,({p-})-E^f (46) 

k 

which only depends on the densities JJ] and their gradients. Inserting this into the density 
equations (p3D finally leads to the reduced equations 



+ ^-4 + ^-^. (47) 



dr 

If we neglect products of spatial derivatives we end up with the coupled Burgers equations 



On the rigth-hand side of this equation we have a sum of diffusion terms with density- 
dependent diffusion functions V^. These cause a smoothing of sudden density changes and 
prevent the formation of shock waves. This is the reason why density gradients and especially 
products of spatial derivatives are normally negligible which justifies the approximation made 
with equation (^8]). Apart from this, the diffusion terms are very helpful for efficient and 
stable numerical integration schemes. 

We will now focus on the simulation of multi-lane traffic. As an example, we investigate 
a two-lane autobahn (i.e. i G {1,2}). For reasons of simplicity, on both lanes the velocity- 
density relations Vf{pi) and pressure relations Vi{pi) will be chosen identically. This is at 
least justified for congested traffic (with a density of 30 vehicles per kilometer and lane 
or more). The corresponding relations are depicted in Figures and ^. They have been 
constructed from empirical data of the Dutch highway A9 between Haarlem and Amsterdam 
with a speed limit of 120km/h and take into account corrections of the traffic equations for 



high densities (for details cf. Refs. [|T^,0). 

The lane-changing rates are chosen in accordance with an empirically validated 



model 50 



1 



± = A^Pi(Pmax - Pi±l) ■ (49) 



Therefore, l/rf" is proportional to the vehicle density pi which refiects the grade of obstruc- 
tion by slower vehicles on lane i. The factor (pmax — Pi±i) refiects that vehicles can change to 



the neighboring lane i ± 1 less frequently the more the density on it reaches the maximum 
density Pmax- For German autobahns the parameters have the following values: 

Pi = 0.176 ■ 10"^ , p2 = 0-056 ■ 10-^ , /3f = = . (50) 

The relation Pi > P2 originates from the fact that the left lane is preferred in Germany, 
since overtaking is forbidden on the right-hand lane. 

Bottleneck situations can be simulated in the following way: We will assume that the right 
lane is closed between places ro and ri. Then, the lane-changing rate Pf will be considerably 
increased, but P2 will be zero on this stretch and already a certain interval Ar before (i.e. 
for ro — Ar < r < ri). P^ and Ar must be chosen sufficiently large so that the right lane is 
empty up to the beginning tq of the bottleneck. Simulation results for the traffic dynamics 
above and below capacity are presented in Figures |^ and ^, respectively. 



V. SUMMARY AND OUTLOOK 

In this paper we have derived a macroscopic traffic model for uni-directional multi-lane 
roads. Our considerations started from plausible assumptions about the behavior of driver- 
vehicle units regarding acceleration, overtaking, deceleration, and lane-changing maneuvers. 
The resulting gas-kinetic traffic model is a generalization of Paveri-Fontana's Boltzmann- 
like traffic equation. It can be extended to situations where different vehicle types or driving 



styles are to be investigated [|T9| 



The gas-kinetic traffic equations not only allow to derive dynamic equations for the vehi- 
cle density on each lane, but also for the average velocity. In this way we were able to extend 
and correct previous phenomenological multi-lane models. Overtaking and lane-changing 
maneuvers are explicitly taken into account, so that the interactions between neighboring 
lanes are included. 

We have then eliminated the velocity equations in order to obtain a reduced model that 
allows efficient computer simulations. The resulting density equations describe the average 
temporal evolution of traffic on a slow time scale. They contain diffusion terms which diverge 
at maximum density Pmax if the finite space requirements of vehicles are taken into account. 



This guarantees that pmax cannot be exceeded and density shocks are smoothed out. The 
latter is important for reahstic results and stable numerical integration schemes. Finally, 
the reduced multi-lane traffic model has been applied to the difficult case of bottleneck 
situations. The computational results were very plausible. Consequently, the model can be 
used to investigate a number of questions concerning the optimization of traffic flow: 

1. In which way does on-ramp traffic influence and destabilize the traffic flow on the other 

lanes? How does the destabilization effect depend on the traffic volume, the length of 
the on-ramp lane, the total lane number, etc.? 

2. In case of a reduction of the number of lanes, is it better to close the left-most or the 
right-most lane? 

3. Is the organization of American freeways or of European autobahns more efficient, or 
is it a suitable mixture of both? Remember that American freeways are characterized 
by uniform speed limits and the fact that overtaking as well as lane changing is allowed 
on both neighboring lanes. In contrast, on European autobahns often no speed limit is 
prescribed (at least in Germany) and average velocity normally increases with growing 
lane number since overtaking is only allowed on the left-hand lane. 

4. In which traffic situations do stay-in-lane recommendations increase the efficiency of 
roads? 
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FIG. 1. The chosen velocity-density relation Vf{pi) ( — ) and the corresponding empirical data 
from the Dutch autobahn A9 (O). 
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FIG. 2. Representation of the density- gradients of the idealized traffic pressure Vi = pOf{p) 
of point-hke vehicles ( — ) and the corrected pressure relation ( — ) which takes into account their 
finite space requirements. Obviously, the increase of the corrected traffic pressure with density and 
the corrected traffic pressure itself diverges at the maximum density /Omax, so that the latter cannot 
be exceeded. For this reason, the diffusion funktions X'fc also diverge for /j^ — > Pmax- The pressure 
relations have been reconstructed from empirical data by means of theoretical relations 10|,19|. 
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FIG. 3. Spatio-temporal evolution of the time-averaged densities on a two-lane freeway stretch 
of 40 km length with open boundary conditions in the case of an overloaded temporary bottleneck 
situation (above: left lane; below: right lane). The right lane is closed between t = lOmin und 
t = 60min on the stretch between ro = 21 km and ri = 25 km, and the vehicles on lane 1 try to 
get on lane 2 beginning at ro — Ar = 20 km. This suddenly increases the density on the left lane in 
the region of the bottleneck, whereas the right lane becomes empty. Already after a short time the 
extremest clustering develops at the beginning of the bottleneck, where the vehicles of the closed 
lane try to squeeze in the left lane. Since the capacity of the remaining lane is smaller than the total 
traffic volume, the left lane becomes overloaded. For this reason a congestion running upstream 
(tailback) builds up on both lanes. On the other hand, the density after the bottleneck, where two 
lanes are available again, is smaller than in front of it so that the vehicles can accelerate there. As 
a consequence, the traffic situation already recovers in the course of the bottleneck. At t = 60 min, 
the lane closure is lifted and the traffic jam disappears. (Note: Due to the different lane-changing 
rates the equilibrium density is somewhat greater than 40 vehicles per kilometer and lane on the 
left lane and somewhat smaller on the right lane.) 
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FIG. 4. Same as Figure |3|, but for a bottleneck situation below capacity. At an average density 
of 10 vehicles per kilometer and lane the traffic capacity of one lane is large enough to cope with 
the total traffic volume. Therefore, no congestion running upstream builds up, but on the left lane 
the density is increased in the region of the bottleneck. This clustering only slowly dissolves in the 
course of the road. After the lane closure is lifted, the traffic jam, which was previously localized at 
the bottleneck, causes a damped density wave. This propagates along the freeway with a velocity 
that is slower than the average vehicle velocity. Due to lane-changing maneuvers, the right lane 
also develops a propagating density wave. 



